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Abstract 

It is well known since the pioneering work of Goldbeter and Koshland [Proc. Natl. Acad. 
Sci. USA, vol. 78, pp. 6840-6844 (1981)] that cellular phosphorylation- dephosphorylation 
cycle (PdPC), catalyzed by kinase and phosphatase under saturated condition with zeroth order 
fvq ' enzyme kinetics, exhibits ultrasensitivity, sharp transition. We analyse the dynamics aspects 

^sO . of the zeroth order PdPC kinetics and show a critical slowdown akin to the phase transition 

ly-j ' in condensed matter physics. We demonstrate that an extremely simple, though somewhat 

mathematically "singular" model is a faithful representation of the ultrasentivity phenomenon. 
' The simplified mathematical model will be valuable, as a component, in developing complex 

cellular signaling netowrk theory as well as having a pedagogic value. 
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1 Introduction 

Cellular informations are stored in the genome; cellular functions, on the other hand, are carried out 
by biochemistry. The communications within and between cells are based on biochemical reactions 
involves enzyme molecules 021. In particular, reversible chemical modification of enzymes, via 
phosphorylation- dephosphorylation cycle (PdPC), is one of the most widely used mechanisms in 
cellular signaling and regulation |@J. In its simplest, canonical form, an enzyme E can be phos- 
phorylated to become E*, and then dephosphorylated to be back to E. These two reactions are 
catalyzed by their respective enzymes called kinase and phosphatase. Biologically, the E form of 
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the enzyme is inactive while the E* form is active. Hence, a change in the concentation of the 
kinase (or of the phosphatase) leads to a change in the activity of the enzyme lfl4l [TTl. 

Mathematical modeling of the PdPC kinetics began in the late 1970s HH. In 1981, Goldbeter 
and Koshland discovered a surprising result: If both kinase and phosphatase, acting as enzymes for 
their respective reactions, are operating under a saturated condition, then the change of E — > E* can 
be extreme sharp in response to a change in the kinase concentration. They called this phenomenon 
zeroth-order ultrasensitivity ||6]. Before that work, biochemists had believed a sharp, sigmoidal 
transition is the signature of only enzyme systems with multiple subunits or monomer with multiple 
binding sites; sigmoidal shaped transition is associated with equilibrium cooperativity. What had 
not been fully appreciated until recent years is that the PdPC is not a closed equilibrium system, but 
in fact an open-chemical system lfl4l . 

There is now a sizable literature on ultrasensitive PdPC and related biochemical reactions. See 
lfl8l [H for two very recent reports. But all the studies, as far as we know, have been limited to 
the steady state analysis, with the exception of |fT9l . The present paper reports a study of the time- 
dependent behavior of ultrasensitive PdPC. We discover that in addition to the sharp transition, the 
kinetic of ultrasensitive PdPC exhibit a critical slowdown, i.e., the time to reach steady state tends 
to infinit when the system approaches the the critical transition point. In a separated study, we 
have also discovered that at the critical point of transition, the ultrasensitive PdPC exhibits large 
fluctuations lfl31[T3"1 . 

Putting all these together, this is a reminiscent of the critical phenonemon in physics ll3l flTTl. 
The classical phase transition problems in condensed matter physics, however, are equilibrium phe- 
nomena. The ultrasensitivity discussed in the present work is a nonequilibrium phenomenon. PdPC 
is a driven chemical reaction system lfl4l . 

2 Phosphorylation-Dephosphorylation Cycle Kinetics 

The canonical kinetic system for the PdPC is (6) |9j 13 

E + K S> EK ^> E* + K, E* +P S E*P ^ E + P, (1) 

fc-ll fe-21 

in which K and P represent kinase and phosphatase. We follow the Michaelis-Menten approach 
to model the enzyme reactions iflOl [TTl. If the amount of kinase is significantlly smaller than that 
of E, and the amount of phosphatase is significantly smaller than that of E* , then according to the 
Michaelis-Menten theory of enzyme kinetics, we have iTTOl ITTl 

dc E * _ Vice V-ic-e* ^ 

dr Kmi + ce Km2 + ce* 1 

in which the total concentration ce+ce* = ct is a constant. Here ce and ce* are the concentrations 
of E and E* respectively. The parameters in Eq. ((2j are related to the rate constants in £T|): 

K Mi = k ~ a + ki2 , (z = 1,2); Vi = k 12 [K] tot , V 2 = k 22 [P] tot , (3) 
kit 

in which [K] to t and [P]tot are the total concentrations of kinase and phosphatase, respectively. 

Equation (0 can be rigorously justified via singular perturbation theory ifTOl fTI. It is valid for 
all times, as long as one can neglects the extremely fast kinetics. Eq. Q in fact is exact for finding 
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the steady state of the reaction system in ©■ The approximation comes from ce + ce* = ct since 
it neglects the concentrations of EK and E*P. 
We introduce nondimensionalized variables 

x=«Z,t = %, (4) 
ct ct 

then we have 

dx 9(1 — x) x 
~dt ~ Ki + l-x~ K 2 + x' 

V\ k\2\K]tot K M \ r , K M 2 ,„ 
U = 77 = 1 — F75i — ' Kl = ' K2 = • ^ 

V 2 k 2 2[P\tot C T C T 

Eq. ([5]) is the starting point of our current study. The steady state x* as a function of 9 and Ks is 
well understood JUITJ. We investigate its time-dependent behavior. 



where 



3 A Simple Mathematical Model for Ultrasensitivity 

The two terms on the right-hand-side of Eq. © represent the Michaelis-Menten rate law of en- 
zyme catalyzed reactions, for the kinase and the phosphatase, respectively. Taking the term r = 
x/(K 2 + x) as example, it has the celebrated "hyperbolic" saturation form, also known as double 
reciprocal since 1 jr is a linear function of 1 jx flTJ . Therefore, the rate r as a function of the substrate 
concentration x is linear, i.e., first-order, when K 2 3> x, but zeroth order when K 2 <C x. When an 
enzyme is operating under zeroth order condition, the catalyzed reaction rate is independent of the 
substrate concentration. 

The important discovery of Goldbeter and Koshland [6] is that when both kinase and phos- 
phatase in (Q]) are operating under zeroth order condition, the steady state concentration of phospho- 
rylated enzyme, i.e., x has a very sharp response to the 6, the activation parameter. 

One can solve the steady state of Eq. (f5]), as a root of a quadratic equation, then let K\ and K 2 — > 
0. This is the standard method of attack. However, if one lets K\ = K 2 = 0, one immediately 
obtains 

1 — X X 

Being mathematically careful with Eq. (0, one has 

x = 0; or x = 1; or 6 = 1, if x ^ 0, x ^ 1. (8) 

Noting that x is a mono tonic increasing function of 9 in Eq. ©, we have steady state x as a function 
of 9 given in Fig. [TJ\. 

This way of setting K\ = K 2 = is in fact the spirit of singular perturbation method iflOl . Fol- 
lowing the same approach and treating K\ = K 2 = e as a small parameter, we have the unperturbed 
equation 

5 = 0-1, (0<z<l) (9) 
at 

The time-dependent solution to the equation is 

x(t) = (9 -l)t + x(0), for 0<x<l. (10) 
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Figure 1: The simple zeroth-order, ultrasensitivity transition in a nutshell: x represents in the frac- 
tion of enzyme being phosphorylated; and 9 is the ratio of the activities of the kinase and the 
phosphatase. (A) The activation curve, steady state x as a function of 9, shows a sharp transition. 
(B) The time-dependent x(t) is a linear function of time t, not exponential as usually expected from 
first-order chemical kinetics. It reaches steady state in finite time too. (Q The time to reach the 
steady state, t^, exhibit a critical slowdown near the midpoint of the transition, when 9 = 1. (D) 
summarizes the key mathematical results associated with Eq. [9] 



This is shown in Fig. Qjj. The most interesting insight from this simple solution is that the time 
reaching the steady state, t^, is finite when 9^1. Furthermore, the time 



g(Q) 

1-6* 

1-1(0) 
6»-l 



< 1, 

(11) 

> 1. 



This is shown in Fig. \Vp. 

The 9 = 1 is known as the mid-point of transition from unphosphorylated E to the phospho- 
rylated E* . It is the "critical point" of the transition. A shape transition like in Fig. [IK is usually 
identified as cooperative phase transition in condensed matter physics. The result on the time in Eq. 
(fTTT) and Fig. \Vp further collaborate this analogue: One of the characteristics of phase transition is 
the "critical slowdown" in its dynamics [17]. In the theory of phase transition, critical slowdown 
and shape transition are also intimately related to large fluctuations ifTTl . This last aspect of PdPC 
has been discussed in Ifl3l[l5l . 

The easiness and transparency of Eq. © makes it very attractive as a simple model for ul- 
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trasensitivity. The remaining of the paper is to provide a more rigorous justification of the results 
obtained from Eq. (©, summarized in Fig. [TjD, in terms of a perturbation expansion approach. More 
specifically, we shall try to provide perturbation corrections to the results in Figs. [T}3 andQJ]. 



4 Exact Solution to PdPC Kinetics Eq. ® 



To justify the above very simple model for ultrasensitivity, one can in fact solve exactly the nonlinear 
differential equation in ([5]). For simplicity, we shall assume K\ = K 2 = K. The Eq. §5$ can be 
re-written as 

dx _ (1 - 9)x 2 + (9 — 9K — K — l)x + OK 
~dt ~ (K + l-x)(K + x) ' 

the numerator of whose right-hand-side has two roots x\ and x 2 : 



(12) 



Xl,2 



-9 + 9K + K+l± y / (l + 2K)(6 - l) 2 + K 2 (9 + l) 2 



(13) 



One of the x\ and x 2 is £ [0, 1], and it is the steady state of the Eq. ((5]>. Let us denote it by x* , then 
the other one is x' = 9K/[x*(l - 9)}. If < 9 < 1, x 2 G [0, 1], and it is the steady state of the Eq. 
©. We also see that K ->• 0, x 2 -> 0. If 9 > 1, x\ € [0, 1], and we have K -> 0, x x -t 1. 

Applying separation of variables and the method of partial fraction to Eq. (fT2l . we have the 
solution to the Eq. ((5]> 



x(0) - x(t) + 



K 



In 



xit) 



x{0) 



X" 



x(t) 



x(0) 



X' 



B 



(14) 



in which, 



A 



1 + K-9K- (l + 9)x* 



D 



1 + K-9K- (l + 0)x' 



From Eq. (fT4l . it is immediately clear that if K = 0, then x = (9 — l)t + x(0). This is 
the result in Eq. (fTOl ). Furthermore, we can compute the "time to steady state". If K — > 0, 



0, = fzf > and for x(t c 



1 ^ _ g(Q)-i 



t — > (xq — x(t)) /(l — 0). Hence for z(ioo) 
These are the results in Eq. (fTTT >. 

For ^ 0, mathematically exponential decay takes infinite time to actually reach the steady 
state. To give a measure of the time, however, we consider the time from x(0) = x* + 0.01 to 
xitoo) = x* + 0.001, for < 9 < 1. Similarly, for 9 > 1, we consider the time from x(0) = 
x* - 0.01 to x(too) = x* - 0.001. 

Fig. |2]shows the as a function of 9, with various K values. It is shown that with sufficiently 
small K, the result in Eq. (|TTT> is indeed valid; There is a critical slowdown at the 9 = 1 for zeroth 
order ultrasensitivity with small K. 

We can also consider the time-dependent solution given in Eq. (fT4l . x(t). Fig. [TJ3 shows simple 
linear functions of time, each of which reaches steady state in finite time. To better understand this 
observation, Fig. [3] shows that the x(t) asymptotically approaches to Eq. (flOT i with diminishing K. 
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0: kinase to phosphatase activity 

Figure 2: Time to reach steady state (x*), defined as from initial value x(0) = x* ±0.01 to x(ioo) = 
x* ± 0.001. ± corresponds to 9 < 1 and > 1, respectively. The values for K are 10 (open circles), 
1 (open squares), 0.1 (filled circles), 0.01 (filled triangle), and 0.001 (filled square). The solid line 
is from Eq. (fTTb - 

5 Discussion 

In developing a mathematical model for a scientific problem, one should always follow Occam's 
Razor which states, in Albert Einstein's words "as simple as possible, but not simpler." The present 
paper contains two results: First, it proposes an extremely simple equation for the zeroth-order 
PdPC, and provides a mathematical proof that the simple model captures all the essential features 
of the well-accepted, but more complex model. Second, we analyzed the time-dependent behavior 
of zeroth-order ultrasensitivity and demonstrated a critical slowdown at the transition point; Thus it 
is a critical phenomenon. 

Critical phenomena are the emergent properties of cooperative processes E [IT] [17). Calling 
the zeroth-order ultrasensitivity a critical phenomenon, thus, begs an answer to the question "what 
is the cooperativity in ultrasensitivity"? Answer to this deeper question is outside the scope of the 
present paper; we refer the readers to two recent papers on the so-called temporal cooperativity 
JEHU. As a matter of fact, a nonequilibrium critical phenomenon in a driven system emerging 
from a cooperativity in time seems to be a novel concept offered by the present mathematical model 
of the biological system. 
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Figure 3: Time dependence of PdPC kinetics, x represents the fraction of substrate enzyme being 
phosphorylated (eg* /or)- All curves have initial value x(0) = 0.6 and 9 = 0.7. Different curves 
correspond to different K values: solid lines from top to bottom correspond to K = 1, 0.1, 0.05, and 
0.01, respectively. The dashed line is from Eq. ([Tol l. With diminishing K, the curves asymptotically 
approaches to the dashed line. 
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